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Abstract 

We initiate a systematic study of the dynamics of multi-particle systems with su- 
persymmetric Van der Waals and electron-monopole type interactions. The static in- 
teraction allows a complex continuum of ground state configurations, while the Lorentz 
interaction tends to counteract this configurational fluidity by magnetic trapping, thus 
producing an exotic low temperature phase of matter aptly named supergoop. Such 
systems arise naturally in N = 2 gauge theories as monopole-dyon mixtures, and in 
string theory as collections of particles or black holes obtained by wrapping D-branes 
on internal space cycles. After discussing the general system and its relation to quiver 
quantum mechanics, we focus on the case of three particles. We give an exhaustive 
enumeration of the classical and quantum ground states of a probe in an arbitrary 
background with two fixed centers. We uncover a hidden conserved charge and show 
that the dynamics of the probe is classically integrable. In contrast, the dynamics of 
one heavy and two light particles moving on a line shows a nontrivial transition to 
chaos, which we exhibit by studying the Poincare sections. Finally we explore the 
complex dynamics of a probe particle in a background with a large number of centers, 
observing hints of ergodicity breaking. We conclude by discussing possible implications 
in a holographic context. 
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1 Introduction 

1.1 String glasses 

In the macroscopic world, where everything flows and nothing abides, mayhem and disor- 
der rule. Most dynamical systems are not integrable but chaotic. Most materials are not 
crystalline but amorphous. They are glasses, eternally in search of their true equilibrium. 
Although glassy systems are still notoriously resistant to fundamental theoretical understand- 
ing, much progress has been made over the past thirty years [3], and remarkable organizing 
principles have been uncovered. A beautiful early example is given by the simplest toy model 
of a spin glass, the Sherrington-Kirkpatrick model, for which the space of low temperature 
configurations organizes itself in a hierarchic treelike fashion [1, 2]. 

String theory presents us with similarly, if not more, complex systems, on microscopic 
scales. These manifest themselves when studying the vast microstates of black holes or the 
vast number of flux compactifications, to give but two examples. Consider for instance (as 
reviewed in detail in [2]), a "3-charge" D4-brane wrapped on a smooth four-cycle £ inside 
a six-torus T 6 = (T 2 )i x (T 2 ) 2 X (T 2 ) 3 , and bound to n pointlike DO-branes. Denote the 
number of intersection points of X with the sub-tori (T 2 )a by P A , A = 1,2,3 — these are 
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the D4-charges of the system. Then E has triple self-intersection product V 3 = 6P l P 2 P 3 , 
and of order V 3 worldvolume deformation and flux degrees of freedom. 

In the regime n>? 3 , the pointlike DO-brane degrees of freedom dominate the degeneracy 
fl of supersymmetric ground states, and one easily computes: 



This agrees with the Bekenstein-Hawking entropy of the D4-D0 black hole. In its 5-dimensional 
uplifted version [4] , the computation is in essence an application of the Cardy formula in the 
context of the AdS 3 /CFT 2 correspondence [5]. 

However, away from the Cardy regime, i.e. when V 3 > n, the order V 3 D4-degrees of 
freedom (deformations of £ and turning on fluxes) become entropically dominant and the 
counting problem becomes far more intricate: Each choice of worldvolume fluxes induces a 
different, highly complex potential on the moduli space of deformations of S, and the ground 
states of the system correspond to the minima of this vast D-brane landscape. On the other 
hand, on the black hole side, an exponentially large number of molecule-like, multicentered 
stationary black hole bound state configurations appears [7], all with the same total charge, 
and they entropically dominate the single centered black hole [8]. In [9], it was argued in a 
probe analysis that this picture persists at least metastably for nonextremal black holes, up 
to a critical temperature where the single-centered D4-D0 black hole regains dominance (see 
also [10] for a closely related analysis). 

This transition between a single black hole and a zoo of metastable multi-black hole con- 
figurations is reminiscent of a structural glass transition and the rugged free energy landscape 
picture associated to the glass phase. Indeed, once in a particular multicentered configura- 
tion, it may take exponentially long for the system to find the most entropic "true equilib- 
rium" configuration, as it proceeds through exponentially suppressed thermal and quantum 
tunneling, and furthermore far from the true equilibrium state, the preferred direction of 
local relaxation processes is likely to push the system into the direction other highly stable 
quasi-equilibrium states instead of the true maximal entropy equilibrium configuration. 

1.2 Supergoop 

When sufficiently far separated and moving slowly close to a ground state configuration, the 
black hole constituents can be thought of as pointlike particles, moving in an approximately 
flat background, interacting with each other through specific static and velocity dependent 
interactions. These effective inter-particle interactions are highly constrained by the fact that 




(1.1) 
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these BPS systems preserve four supercharges: A nonrenormalization theorem [6] implies 
that once a metric has been fixed on the configuration space, the static and first order 
velocity dependent interactions are of a fixed form. For the flat metric: 



^— ' 2m. 

p=i 

where A p is the vector potential produced at x p by a collection of Dirac monopoles of charge 
n pq at positions x g . The coefficients K pq = —K qp equal the electric-magnetic symplectic 
products between the charges of the centers p and q. The parameters 9 P and masses m p are 
fixed by the BPS central charge of each center. If the configuration space metric is not flat, 
the m p may depend on x p . 

As a result of this nonrenormalization theorem, exactly the same supersymmetric La- 
grangians also appear in very different physical contexts where four supercharges are pre- 
served and the low energy degrees of freedom can be identified with spatial positions. One 
example is a mixture of well-separated elementary particles obtained by wrapping D-branes 
around various internal cycles of a Calabi-Yau manifold, interacting with each other through 
gravitational, vector and scalar interactions. Clearly this can be viewed as an extreme limit 
of the multi-black hole systems considered above, where the dyonic black holes have been re- 
placed by dyonic particles. Another example are monopoles and dyons in M = 2 Yang-Mills 
theories [15]. 

A more remote example [6] is a collection of space-localized wrapped D-branes at weak 
string coupling in the substringy distance regime. Their low energy degrees of freedom are 
given by a (0 + 1) -dimensional supersymmetric quiver quantum mechanics, with a position 
3-vector and a U(l) gauge symmetry for each singly wrapped brane (identified with the 
nodes of the quiver) and the lightest brane-brane stretched open string modes represented 
as bifundamental oscillator degrees of freedom (identified with the arrows of the quiver). 
When the branes are all well separated, i.e. when the quiver theory is on the Coulomb 
branch, the open string modes become very massive and can be integrated out. Again, the 
resulting effective theory for the position degrees of freedom must necessarily be of the form 
(1.2) fixed by supersymmetry. The coefficients n pq are now identified with the net number 
of arrows between two nodes. 

Thus this type of supersymmetric multi-particle mechanics appears in many contexts, 
in widely different regimes. Much effort has been put into understanding and counting the 
supersymmetric ground states of such systems, in part because of their key role in physics 
derivations of BPS wall-crossing formulae [6, 8, 16, 17, 18]. However, little has been said 
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about excitations or dynamics for these systems. There are a few exceptions: [19] studied 
the classical and quantum dynamics of the two-particle system and found it was integrable, 
and in [9] the persistence of the black hole molecular configurations at finite temperature 
was studied. However, no studies of multi-particle dynamics or statistical mechanics have 
been done so far. In this paper we wish to take steps in these directions. 

Besides the motivation for understanding D-brane and black hole statistical mechanics 
and their potentially interesting interpretation as holographic glasses, such studies would also 
be of intrinsic interest, as these systems are rather unusual in several aspects. Due to the 
special form of the potential (1.2), an iV-particle bound state will have a 2(iV— l)-dimensional 
moduli space of zero energy ground state configurations folded in a very complicated way 
into the 3(iV — l)-dimensional full configuration space (factoring out the center of mass). 
Naively one might therefore think that even at very low temperatures, the system would 
behave like a liquid, exploring large parts of the configuration space by flowing along the 
continuous valley of minimal energy configurations. As a simple example consider the case 
N = 2. The particle distance is fixed and the moduli space is a sphere. One might think a 
probability density initially localized near a point on that sphere would quickly diffuse out 
over it. However, due to the effective electron-monopole Lorentz interaction between the 
particles, this is not quite right, as diffusion is obstructed by magnetic trapping. Another 
way of understanding this is conservation of angular momentum: Monopole-electron pairs 
carry intrinsic spin directed along their connecting axis, of magnitude equal to half their 
symplectic product. Hence they behave like gyroscopes. They resist changing direction; 
kicks will just cause them to wobble. 

Thus it is natural to hypothesize that these supersymmetric multi-particle systems behave 
partly like a liquid and partly like a solid at very low temperatures — like goop. We will 
therefore refer to this peculiar state of matter as supergoop. 

1.3 Dynamics 

Many times we study Hamiltonian systems that are classically integrable. For this to be 
the case one requires the existence of at least iV conserved charges (with all mutual Poisson 
brackets vanishing) for a system with a 2A-dimensional phase space. Examples include single 
one-dimensional particles with an arbitrary potential, since the energy is conserved, and two 
body problems with a central force. Phase space trajectories of a classically integrable system 
will map iV-dimensional tori. Generally, however, our system will not be classically integrable 
and we must confront a chaotic system. The simplest example of a chaotic system is the 
double pendulum which has a four-dimensional phase space with a single conserved quantity: 
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the energy. One then studies the phase space trajectories of the double pendulum as a 
function of increasing energy. For sufficiently low energies the trajectories are constrained 
to live on a two-dimensional torus displaying quasi-integrable behavior. As the energy is 
increased this torus is deformed and eventually breaks apart into smaller tori. This process 
is seen to continue until there is no visible structure in the phase diagram, i.e. the system 
tends toward ergodicity. What is perhaps most remarkable about the transition to chaos is 
that it occurs in a gradual fashion in which smaller and smaller islands of regular behavior 
are spawned before the system loses all manifest structure. A crucial question, especially for 
a system with many degrees of freedom, is quantifying when all ordered behavior disappears 
and how it depends on the parameters of the system (see for example [20, 21]). 

It is our aim to begin a systematic study of the dynamical aspects of the underlying 
brane system on the Coulomb branch. In this paper we mainly address the case of a three 
particle configuration. This is already a difficult non-integrable three-body problem. To 
render the problem tractable, we study first the classical ground states and subsequently the 
motion of a probe particle in a fixed background consisting of a two-centered bound state. 
Remarkably, we discover that the motion of the probe is classically integrable! This is due to 
the presence of an additional hidden conserved quantity and is somewhat reminiscent of the 
integrability of a Newtonian probe particle interacting gravitationally with a background of 
two fixed masses, as discovered by Euler and Jacobi. We then study the transition to chaos 
for a system of two probes in the presence of a heavy fixed particle, with the motion re- 
stricted to live on a line. This setup is directly analogous to the double pendulum, allowing 
us to exploit many of the tools developed for the study of the double pendulum. As for 
the double pendulum, by studying Poincare sections we observe the formation of islands in 
phase space and the eventual transition to global chaos with no apparent structure in phase 
space. Finally, we begin to address the far more intricate dynamics of a system with a large 
number of centers. We provide a brief exposition of the behavior of a probe particle inside 
a molecule with a given number of fixed centers. We observe highly complex trajectories 
that become trapped for long times in the sense that they do not explore the entire molecule. 

Note added: After submission, we became aware of [11, 12, 13, 14], which elegantly prove 
the classical integrability that we discuss in Section 5. 

2 General Framework 

Consider a system of branes wrapped on the cycles of a six-dimensional compact space, such 
that they are pointlike in the non-compact (3 + 1) -dimensions. The interactions between 
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them are governed by strings whose ends reside on the branes themselves. The low energy 
physics is governed by an M = 4 supersymmetric quiver quantum mechanics [22] . The nodes 
of the quivers have gauge groups associated to them and the low energy string degrees of 
freedom are chiral multiplets transforming in the bifundamental between two given nodes. 
The position degrees of freedom x p of the branes are scalars in the vector multiplets of the 
gauge groups. It was shown in [6] that if we study the Coulomb branch of the branes, i.e. 
integrate out the massive chiral multiplets, a non-trivial potential is generated which governs 
the dynamics of the x p . The system of branes allows for a large family of bound states with 
fixed equilibrium distances as classical ground states. If the branes move too close the 
stretched strings become light and even tachyonic and the system enters the Higgs phase. 1 
On the other hand, if there is a sufficiently large number n# of branes placed at a single point, 
such that the product of the string coupling g s and % becomes large, the system is best 
described by closed string exchange and hence supergravity. In what follows we will consider 
the theory of supersymmetric multiparticle mechanics describing the Coulomb branch of the 
brane system. This has the advantage of being a simple setup which is interesting in and 
of its own right while reproducing many of the features of the multicentered configurations 
that exist in supergravity. 

2.1 Supersymmetric multiparticles 

The theory we consider is the multiparticle supersymmetric mechanics, which we refer to as 
supergoop, studied for example in [6, 7, 19, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32]. We simply 
state the Lagrangian of the system, referring to [6] for details. 
The supergoop Lagrangian is given by: 



where X p is the fermionic superpartner to x p . The D p fields are auxiliary non-dynamical 
scalars. We have defined the functions: 



± 2 p + D 2 p + 2i\ p \) + (~ U p D p + A p ■ *p) + ( C pi X p X i + C pi ■ \°\) > 
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In fact, the weak string coupling limit g s — > always pushes the system to the Higgs phase. 
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is the vector potential of a single magnetic monopole of unit charge at the origin and r pq = 
x p — x g . For the Lagrangian to be supersymmetric we further require K pq = —n qp . In the 
quiver quantum mechanics context, the the number of bifundamentals connecting 

two nodes. The supercharges are given by: 

Q a = - * 'ZX, ■ (Pp - Ap) + KU P , Q P = Y, 1 < X l ■ (P* - A,) " ~KU q . (2.4) 
p q 

The Weyl spinors obey (Ac)* = A a and the cr^ are the usual Pauli matrices. The Hamiltonian 
of our system is defined as: 

H = ^2p p -± p - L , p p = m p ± p + A p . (2.5) 
p 

Upon integrating out the the Z)-terms we find: 

H = £ Vpp ~ Ap) 2 + u p] + E |5r r - ■ V CT V , (2.6) 

P p p<q PI 

where X pq = X p — X q . Notice that the particle interactions include velocity dependent forces. 
Furthermore, the system has three-body interactions due to the appearance of U p in H. 

As usual our Hamiltonian H is related to the supercharges as {Q a , Q i3 }d.b. = —2i5^H. 
However this is most easily checked in the quantum mechanics context where p p — > —iV p 
and we replace the above Dirac bracket relation with the anticommutation relation: 

{Q a ,Q?} = 25iH , (2.7) 

where {A£, A^} = m; 1 ^. 

2.2 Classical features and multicentered black holes 

To study the classical properties of this theory, we can turn the fermionic fields off and only 
consider the bosonic part of (2.1). Static BPS configurations occur when U p = for all p, 
i.e. when: 

A i pa 

Taking the sum over p of (2.8) we find that the 9 p s must satisfy: Yl p @p = 0- The solutions to 
equation (2.8) are bound states of particles, see for example figure 2.1. Given that a system 
of iV particles has 3iV-degrees of freedom which are constrained only by (N — 1) equations, 
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Figure 2.1: Examples of ground states for 100 electric T e = (0, 1) plus 100 magnetic T m = 
(1, 0) particles. 

the bound states have a (2N + 1) -dimensional classical moduli space M.. This moduli space 
cannot be accessed dynamically at low temperatures however, due to the velocity dependent 
forces which constrain the particles to oscillate about a fixed location if given a small kick 
(just like an electron in the presence of a magnetic field). As the energy is increased, the 
rigid structure of the bound state is deteriorated and eventually lost completely. In the 
following sections we study the dynamical features of this system, with a particular focus on 
the three-particle system. 

Equation (2.8) is a familiar expression in supergravity. It recreates the integrability 
condition of [7] for multi-centered black hole bound states in four- dimensional Af = 2 super- 
gravity: 

The above expression involves electric-magnetic charge vectors T = (P 1 , Qi) with a duality 
invariant symplectic product given by: 

(r,f) = P I Q I -Q I P I . (2.10) 

Expression (2.9) also involves a function called the central charge Z p (z) which depends on 
the vector multiplet scalars z a of the supergravity theory, as well as the charge vector T p . 
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At spatial infinity we can write Z p \ r=OQ = m p e iap where m p is the ADM mass of a BPS 
particle of charge T p . If we denote Z = Yl p Z p then the parameter a in (2.9) is given by 
a = arg (Z\ r=00 ). Thus we may rewrite (2.9) as 

J2 ffzi£?> = mp sin ( Qp _ a ) _ (2.11) 

q ^ r PQ 

The supersymmetric multi-particle mechanics in (2.1) may be considered as a toy model for 
the dynamics of the multi-centered black hole bound states if we make the identifications 

K pq = (T p ,Tg) and 9 P = — m p sin(a p — a) . (2-12) 



2.3 Three particles 

Much of the discussion that follows will concern the three particle system, which already 
exhibits rich dynamic and non-dynamic features. Here we describe some of the characteristic 
features of its zero energy configurations. 

Classically, the (supersymmetric) ground states are found by setting U\ — U2 — Us — 0. 
Explicitly: 

Kl2 K31 -e 1 , -^ + ^ = -e 2 , -^ + ^ = e 3 , (2.13) 



2r i2 2r 3 i 2r 12 2r 23 " 2r 31 2r 23 

with # 3 = — (61 + 62). Notice that the third equation follows from the other two. We also 
require that the three relative distances 7*12, r 23 and ri 3 satisfy the triangle inequality. The 
three particles have nine position degrees of freedom and the above equations only constrain 
two of them. Factoring out the center of mass leaves us with (9 — 2 — 3) = 4 unconstrained 
degrees of freedom. Hence, even when three-particle bound states form, there is an infinite 
classical moduli space of connected (and possibly also disconnected) ground states. In the 
case of a two-particle bound state the classical moduli space is simply a two-sphere of fixed 
radius K\z/29i. 

If we have k%\, k 12 and k 23 positive we find there exist scaling solutions given by — > Xk^ 
with A —7- with the k's obeying the triangle inequality [8]. Thus, in this regime the particles 
can come arbitrarily close to each other with no cost in energy. 2 Away from the scaling regime 
the solution to (2.13) corresponds to a bound state for which the particles may oscillate about 



2 When thinking about scaling solutions in the gravitational context from the point of view of a far away 
observer, the scaling solutions are continuously connected to a single center and will look like a single centered 
black hole or particle. On the other hand, nearby observers will observe that the proper distance between 
the particles never shrinks to zero due to the formation of infinite throats (at least at zero temperature). 
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a fixed equilibrium radius upon small perturbations. 

As an example, when 6*3 = —3, 62 = —2, Q\ = 5, K12 = — 1, K43 = — 1 and K23 = 1, 
with the three particles living on a line and particle 3 between particles 2 and 1, we find the 
solution: 

ri2 = 55 (7 + v 7 ^) « 0.57, r 13 = ^ (8 - v 7 ^) ~ 0.12 , r 23 = ^ (l + v^) « 0.45 . 

(2.14) 

Note that the triangle inequality is saturated since we have considered a collinear example. 
We should note, however, that there are clear instances where no solutions exist, such as 
when (K13, k 2 3, #3) > (0,0,0), for example. 

2.4 Regime of validity 

Since we are free to choose the set of a p , (2.12) does not really constrain the values of the 8 P 
in any way. There is, however, a restriction stemming from the requirement of the validity 
of the Coulomb branch description assuming our system comes out of integrating strings 
[6, 34]. The distances between particles must be smaller than the string scale, but larger 
than the ten dimensional Planck scale. For larger distances, the suitable description is given 
by the exchange of light closed strings (in which case supergravity is the reliable description) . 
Furthermore, the velocities should be small compared to the speed of light to avoid higher 
derivative corrections to the non-relativistic Lagrangian (2.1). 

Consider first two particles with masses mi and m 2 and call the string length l s . The 
restriction is found to be [6]: Ki 2 <C l s 9i and \i l 2 s 0\/ ' k\ 2 , where \i = mim 2 / (mi+m 2 ) is the 
reduced mass of two particles. Indeed, if K\ 2 3> l s Q\ the distance between the two particles 
in a bound state will be much larger than l s and the appropriate description becomes that 
of supergravity. When \i ^i/ K i2) °P en strings between the branes become light and 
the appropriate description becomes that of the Higgs branch and eventually the fused D- 
brane system itself. For the multiparticle system, the Coulomb branch description is reliable 
so long as the inter-particle distances are sufficiently large that the massive strings can be 
reliably integrated out, i.e. ^> l s y/\6jjrrii — 9jjrrij\ and sufficiently small that we remain 
in the substringy regime, i.e. l s . 

3 Classical Phase Space 

Having discussed the general framework of the system under study, we now discuss some of 
its dynamical features, beginning with the classical phase space. Recall that the Hamiltonian 



11 



of our system is given by: 



p v 

The Hamilton equations of motion are given by: 

\7 Xp H = -p p , \7 Pp H = ± p . (3.2) 

For iV-particles we have a 3 x 2N = QN dimensional phase space. As manifest conserved 
quantities we have the energy, the center of mass momentum, and the center of mass angular 
momentum. 

3.1 Two particles are integrable 

We review the classical properties of the two particle problem in appendix A. Recall that a 
classically integrable system with a 2iV-dimensional phase space has at least N conserved 
quantities with mutually vanishing Poisson brackets. Phase space trajectories for integrable 
systems reside on N- dimensional tori. In the case of two-particles we have a 12-dimensional 
phase space. There are six manifest conserved quantities given by the net momentum and 
angular momentum. As shown by D'Hoker and Vinet [25, 26], the presence of a conserved 
Runge-Lenz vector (A. 2) leads to an enhanced 5*0(3, 1) symmetry. The angular momentum 
and Runge-Lenz vector are three-vectors, the Hamiltonian is a scalar and there exist two 
relations amongst the seven quantities, hence there is a total of (3 + 3 + 2) = 8 conserved 
quantities. Factoring out the center of mass yields a (maximally) super-integrable system. 3 
The super-integrability implies its equations can be separated in more than one coordinate 
system and one can solve for the quantum mechanical spectrum algebraically, as done in 
[19]. Further, it implies that trajectories in coordinate space follow paths which are closed 
in the case of bound orbits. 

3.2 Three particles are chaotic 

In the case of three-particles we have an 18- dimensional phase space and there is no longer 
a sufficient number of conserved quantities to render the system integrable. Thus, such 
systems will exhibit chaotic behavior. 

We may use several numerical tools to analyze the chaotic nature of such a system. For 

3 A super-integrable system [35] is a system with a 2iV-dimensional phase space which has more than N 
conserved quantities. A maximally super-integrable system has 2N — 1 conserved quantities. 
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instance, we can study the Lyapunov exponent A parameterizing the divergence of phase 
space trajectories with nearby initial conditions. Given two trajectories in phase space with 
initial separation 5zq, the Lyapunov exponent is defined by the limit: 

A = lim lim - log ^ . (3-3) 

t^oo <5zo->0 t OZq 

We can also study Poincare sections in phase space. These are found by recording the 
location of a trajectory in a particular subspace of phase space each time it crosses some 
fiducial point (such as crossing the origin with positive velocity). These are particularly 
useful for lower dimensional systems such as the double pendulum, where they clearly depict 
the breakdown of the integrable motion on a two-torus as the energy is increased (see chapter 
11 of [20] for a discussion). We will discuss and examine the Poincare sections of a collinear 
three particle system in section 6. 

Six of the phase space dimensions can be eliminated from net momentum conservation 
and factoring out the center of mass. We can also kill another four due to net angular 
momentum and energy conservation. The remaining 8-dimensions in phase space (as far 
as we know) are unconstrained by symmetries. Needless to say, systems with more than 
three particles will also display chaotic properties. A simpler setup, which we refer to as the 
Euler-Jacobi setup, is that of a probe particle moving in the background of two fixed centers. 

4 Euler-Jacobi Ground States 

The simplest question we can ask about our system is what the (supersymmetric) ground 
states are, both classically and quantum mechanically. Classically there may be continuous 
moduli spaces of zero energy configurations. Quantum mechanically, given that the probe is 
a particle in the presence of a background magnetic field, we expect the continuous classical 
moduli space to give rise to a degenerate set of quantum ground states due to Landau 
degeneracies. 

4.1 Euler-Jacobi three body problem 

We will consider a probe particle of mass m 3 in the background of two fixed centers unless 
otherwise specified. The background particles have masses m\ and rri2 both very large 
compared to 777.3, charge vectors Ti and T 2 with symplectic product K12 > and Fayet- 
Iliopoulos constant Q\ = —6 2 < 0. They sit along the z-axis at z = ±Ki 2 /4#i = ±a. By 
choosing mi and 7772 much larger than 777,3, the backreaction of the probe on the fixed centers is 
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suppressed by 0(m 3 /mi,m 3 /m 2 ). The probe also has charge vector 73 and Fayet-Iliopoulos 
constant 83. 

We may also consider the possibility of forming supersymmetric bound states between 
the probe and the fixed centers. In such a case, a non-zero 63 requires us to modify the 
background condition 8\ = —6 2 , since now the #'s must add to zero. We thus demand 
I $3 1 *C \8± |, |^2| such that the correction to the positions of the original fixed centers is of 
order 0(83/61,83/62). The intersection products of the probe with the centers are given by 
K31 and K32. To avoid any large backreaction from the probe due to the k interactions we 
further require that K31/V31 and K 32 /r 32 are small compared to k 12/7*12 ~ # 2 . 

The Hamiltonian governing the dynamics of the probe is given by: 

JWe = ^" (P3 " A3) 2 + ^- U + ^ + ^) 2 ■ (4.1) 

2m 3 2m 3 \ 2r 31 2r 32 / 

Notice that the two scaling transformations: 

(p 3 (t), x 3 (t) ; K i:j , 9 3 , m 3 , r 12) t) -»■ (<r p 3 (t), A o^ 1 x 3 (t) ; A a 6» 3 , a 2 m 3 , A a -1 r 12 , A t) , 

(4.2) 

generate a family of solutions parameterized by A and a . Given a solution to the equation of 
motion for some 0(1) parameters, we can exploit the scaling symmetries to map the solution 
to a rescaled one in the regime of validity for the Coulomb branch description as discussed in 
section 2.4. In particular, we require A ^> \fa and A o which can be achieved for large a. 
Notice that in this regime the velocity (which scales as cr _1 ) becomes parametrically small. 



4.2 Classical Ground States 

As noted in (2.8), the space of classical ground states M. is given by setting U p = 0. Satisfying 
this condition gives rise to time independent classical bound states. In the probe limit, where 
the two background centers are fixed, this amounts to solving the algebraic equation: 

1 K31 ' 1 K32 , (4.3) 



2 VP 2 + (2 " a) 2 2 + ( z + a) 2 

where p 2 = x 2 + y 2 and <ft = tan _1 (y/a;). One can easily prove that the effective magnetic 
field B = V 3 x A3 is always perpendicular to the tangent of Ai. This remains true for the 
moduli space of a probe in a background of more than two centers as well. 
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Consider first the case 63 = 0. We find: 

p 2 = (a-^) 2 ^ 2 - {a + z) 2 K 2 31 ^ ^ 

( K 31 ~ ^32) 

For p(z) above to have solutions we choose K31 > > K32 and furthermore |«3i| > | 1 , we 
find a continuum of solutions between z — [z-, z+] where: 



z± = ±a 



|«32| 


=F 


«31 


«31 


± 


|«32| 



(4.5) 



Note that z_ < —a < z + < a and thus the #3 = surfaces enclose the fixed charge at z = —a 
in this case. Similarly, for |«3i| < | 1 the probe encloses only the center at z — a. For 
^31 = — ^32 and 83, the classical moduli space M. becomes the z = plane. 

For 6*3 7^ finding p(z) amounts to solving a quartic equation. In order to do so, it is 
convenient to go to prolate spheroidal coordinates: 

p = a^/(e-l)(l-V 2 ) , z = a(iv, = 0, (4.6) 

such that: 

2a9 3 (r] 2 - £ 2 ) = («ai + K 32 )e + (« 31 - ^2)V • (4.7) 

In the above we have implicitly used that 77 G [—1, 1] and £ e [1, 00]. We can easily find a 
solution for 77 = r/(£): 

V(0 = \{S 1 -S 2 ± y /(5 1 -5 2 y + 8(5 1 + 5 2 ^+lQe) , (4.8) 

where Si = K31/ (a9 3 ) and 5 2 = ^32/(0^3)- In figure 4.1 we show the different qualitative 
types of Ai as a function of Si and S 2 . The qualitative features of each region are shown in 
figure 4.2. Notice that upon defining the prolate coordinates (4.6) we have scaled out the 
distance r 12 = a between the fixed centers. To obtain physical distances we simply multiply 
by r 12 = a. 



4.3 Quantum Ground States 

From the classical point of view, our particle is nothing more than a charged particle in the 
presence of magnetic fields constrained to live on a surface. Thus, given a time independent 
supersymmetric bound state configuration we can compute the lowest Landau degeneracies 
d,L by computing the degeneracy of states with vanishing energy for the constrained particle. 
Such a setup has been addressed for non-uniform magnetic fields everywhere normal to the 
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6z 6i 

Figure 4.1: Left: Classical moduli space Ai in the 8\ — <5 2 plane for #3 7^ 0. The nature of 
Ai for the different regions is shown in figure 4.2. Right: Classical moduli space for Ai with 
#3 = 0. In regions i and ii the centers at z = a and z = —a are enclosed respectively. 

surface [36], which is precisely the situation we find ourselves in. Following [6], the lowest 
Landau degeneracies are given by the total flux through the classical moduli space Ai. For 
instance, as we discuss below, the Landau degeneracy of the fixed background is given by 
K2i- Upon studying the phase diagram and corresponding Ai in figures 4.1 and 4.2, we find 
that the total degeneracy is: 

I, II and III : d tot = k 12 x |« 31 | or d tot = k 12 x |« 32 | , (4.9) 
IV and V : d to t = «12 x |« 3 i + k 32 \ ■ (4-10) 

For regions I, II and III, the degeneracy of states depends on which of the two background 
centers is encircled by Ai. Notice there is a jump in the number of ground states as we vary 
5i and 82 ■ Since we are in the probe limit, we expect these results to be correct up to order 
0(«3i/«ia) and 0(k 32 /k 12 ). 

From the supersymmetric quantum mechanics point of view, recall that {Q a , Q^} = 
25^H. In the absence of the probe, the ground state of the background is given by [6]: 

|b> = ^f a (x x - x 2 )l a \0) , A = A x - A 2 . (4.11) 

The center of mass coordinate xq = (miXi + m 2 x 2 ) /{mi + m 2 ) and center of mass spinor 
A = (777.1A1 + m 2 A 2 ) /(mi + m 2 ) drop out and thus |b) is naturally a function of the relative 
background position vector and spinor. The state |0) is annihilated by A and defines a 
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Figure 4.2: Classical moduli space in the 5i — 62 plane for # 3 7^ 0. The order of the figures 
left to right starting at the top are the regions in figure 4.1. 
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Figure 4.3: Three node quiver with a closed loop (left) and without a closed loop (right). 

three-dimensional Hilbert space through action of A. There are K42 ground states filling a 
spin-(/«i2 — l)/2 multiplet. From the last term in the Hamiltonian (2.6), we observe that 
there exist spin-spin couplings between the background spinors Ai and A 2 and the probe 
spinor A3. It is convenient to introduce the relative spinors A 13 = Ai — A 3 and A23 = A 2 — A 3 
and their associated vacua IO23) and IO13), such that A23IO23) = and so on. Both IO23) and 
I O13) have an associated three-dimensional Hilbert space and the general state must be a 
tensor product of all linear combinations of all such states, finally tensored with |b). Given 
that the probe is sensitive to the background B-field, it will go into spin one-half states of 
I O13) and IO23) aligning with the B-fields from the fixed particles at z = ±a. This will split 
the lowest Landau degeneracies. It would be interesting to compute the explicit ground state 
wavefunctions. 

For the sake of completeness we briefly mention another method to compute the number 
of ground states. One can associate a quiver diagram Q to the data (k^, 9j) of a particular 
configuration [22, 6]. It turns out that the dimension of the moduli space M.(Q, N, 9) of 
the quiver Q can be related to the number of BPS ground states. In particular, for the 
three body problem where each particle is a different species we have a quiver theory Q 
with N = (1,1,1), K12 arrows between nodes 1 and 2, « 13 arrows between nodes 1 and 3 
and K23 arrows between nodes 2 and 3. The quiver diagram is presented in figure 4.3. The 
Fayet-Iliopoulos constants 9 V are additional parameters associated with each node v. Ground 
state degeneracies for similar setups to the one we are studying have been computed in [8]. 
Notice that scaling solutions can occur only for quivers with closed loops. In our problem, 
with K12 > 0, we find that regions I, II and III correspond to quivers with closed loops and 
regions IV and V correspond to quivers with no closed loops. 



18 



5 Euler-Jacobi Dynamics: classical integrability 



The equations governing the probe are dictated by the Hamiltonian in (2.5). There are two 
obvious constants of motion in this problem, namely the energy and the angular momentum 
in the direction of the line where the two centers are placed. If the system is to be rendered 
integrable, there must exist a third constant of motion. Such a constant of motion was found 
for the problem of a Newtonian probe interacting gravitationally with a background of two 
fixed massive particles [37, 38], also known as the Euler-Jacobi three-body problem. We will 
show that the analogous problem in the theory under consideration is also integrable. This 
was previously shown and discussed in [11, 12, 13, 14]. 



5.1 Setup and coordinate systems 

Recall that we are considering two fixed background centers sitting on the z-axis at z — 
±K 12 /4a6' 1 = ±a. Let us go to a cylindrical system with metric: 

ds 2 = dp 2 + p 2 d(j) 2 + dz 2 , (5.1) 

where the Cartesian and cylindrical coordinates are related by x = pcos0, y = psin0 
and z — z. One observes that the Lagrangian and Hamiltonian are independent of the 
coordinate which implies a symmetry. The conserved quantity of this symmetry is given 
by the angular momentum in the z-direction, such that the canonical momentum p^ = I is 
constant. The probe Hamiltonian (4.1) in cylindrical coordinates becomes: 

H probe = (p5 - Hp, 0, *)) 9 ij & - Mp, <t>, *)) + {U i P m f 2 > ( 5 - 2 ) 

where the p' { are the conjugate momenta in the cylindrical coordinates. The relation be- 
tween conjugate momenta between the primed and unprimed coordinate systems is pi = 
p'j dx' j /dx l . 

The third constant of motion is not manifest in cylindrical coordinates. One must go to 
the prolate spheroidal coordinates (4.6) with metric: 



ds 2 =a 2^2_ v 



d£ 2 drj 2 



+ 



+ a 2 (e-l)(l-V 2 W ■ (5.3) 



_(e-i) (l-V 2 ). 

Once in this coordinate system, we note that our Hamiltonian takes the following form: 

£2 _ ^2 
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Figure 5.1: The Euler-Jacobi flower. The red balls represent the fixed background centers 
and the blue line represents the classical trajectory of the probe. In this case, the trajectory 
precesses around only one of the fixed centers. 

where depends only on £ and p| and H v depends only on 77 and p 2 . Thus, we can write: 

Hprobe^ 2 — H$ = H v + H probe rj 2 = G , (5.5) 
where G must be a constant of motion. More explicitly we have: 

H = 2 (£ 2 ~ 1) Pi £(«3i - « 32 ) (K31 ~ ^2? 3 £(k 3 i + k 32 + afl 3 fl 

5 Pf 2a 2 m 3 2a 2 m 3 (£ 2 - 1) ^2a 2 m 3 (£ 2 - 1) 8a 2 m 3 (£ 2 - 1) 2am 3 

and 

H = 2 ^ -V 2 ) ! P\ ! v( K 3i + ^32) | («3i + /t 32 ) 2 | g 3 7y(K 3 i - k 32 - a03V) 

v ^ v 2a 2 m 3 2a 2 m 3 (l — rf) ^^2a 2 m 3 {rf — 1) 8a 2 m 3 (l — ?] 2 ) 2am 3 

We conclude that the probe-two-center problem of supergoop is integrable, providing another 
example to the distinguished list of integrable classical systems. In this system, one observes 
highly symmetric spatial trajectories, as illustrated in figure 5.1. 



6 Beyond Euler-Jacobi: the stringy double pendulum 

If we move away from the probe approximation and allow backreaction with the fixed centers, 
our system is no longer integrable and begins to show chaotic features. For instance, one can 
study trajectories in phase space and see whether they are closed. One could also compute 
the Lyapunov coefficient of the system. In figures 6.1 and 6.2 we demonstrate the trajectories 
in phase space for the probe orbiting around both centers as we exit the probe limit. 

As we increase the number of degrees of freedom, the analysis of chaotic systems becomes 
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Figure 6.1: Examples of closed phase space trajectories in the integrable probe regime. The 
plots show slices of phase space in the Cartesian coordinate system. 




increasingly challenging. The canonical example of a double pendulum, which already dis- 
plays a significant set of features generic to chaotic systems at large, can be effectively 
analyzed with the use of Poincare sections. For a double pendulum, the phase space is four- 
dimensional with a single constant of motion - the total energy. Hence the system is not 
integrable. 

6.1 Collinear dynamics 

Away from the probe limit, as we already noted, our system has an 18-dimensional phase 
space and becomes a complicated three-body problem. In order to study the transition to 
chaos of our system, it is instructive to find a setup that allows us to use the same tools 
used to analyze the double pendulum. This can be achieved by restricting the particles to be 
collinear, i.e. placing them on a line and only considering deformations along this direction. 
Notice that a system consisting of particles on a line will stay on the line so long as the 
velocities of the particles are parallel to the line itself. This is because the magnetic force 
v x B will vanish in this situation. Hence, the collinear system is a consistent truncation of 
our original Lagrangian (2.1). This is no longer true for the coplanar case. 

As was already discussed, we need at least three particles to find chaotic features. Three 
backreacting particles on a line have six degrees of freedom with a conserved energy, a 
situation closer to the triple pendulum. We can however take the mass of one of them to be 
much larger than the other two such that they behave as two probes in a fixed background. 
The probes are allowed to interact with each other since we do not restrict the 9i and Kij in 
any way, except £\ 9{ = 0. The equations governing the two probes can be extracted from 
the two-probe Hamiltonian: 

= A + _?L + J_ (e 2 + J* + ««-) 2 + J- (k + J^ + ^lY. (6.i) 

2m 2 2m 3 2m 2 \ 2x 21 2x 23 J 2m 3 \ 2x 31 2x 32 J 

The above Hamiltonian is a good approximation for the three-particle system in the limit 
where m\ 3> m 2 ,m 3 . In this limit one particle becomes non-dynamical and the energy is 
fully contained in the motion of the two light particles. Thus there is a conserved quantity 
associated to the motion of the light particles and we are left with a three-dimensional phase 
space, which is also the dimensionality of the double pendulum phase space. 

The ground state (^ij^jsi) is found by setting U 2 = U 3 = 0. In addition to imposing 
the triangle inequality to fully specify the ground state, we must also declare the ordering 
of the three particles on the line. A given ground state is mapped to a family of ground 
states via the scaling relation (x^x^; k^) — > A^ij^jjii Slightly increasing the energy 
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leads to small oscillations about the equilibrium position {x* 21 ,x% il ). The linearized normal 
frequencies are the eigenvalues of the Q 2 matrix: 



The derivatives of H co i are evaluated at the equilibrium point. Though the general formulae 
for the normal frequencies are quite involved they are readily computed. As an example, the 
normal frequencies for k 13 = k 12 = — k 23 = 1, d\ = — 9 2 = — 1, #3 = 0, and m 2 = m 3 = 1 are: 



6.2 Poincare Sections 

When the system is integrable or quasi-integrable, i.e. for sufficiently low energies, the 
trajectories in the four-dimensional phase space will reside on two-dimensional tori, since 
the system is simply given by two linearly coupled oscillators. Since it is hard to visualize 
motion on the torus, we study instead particular snapshots of the system, known as Poincare 
sections (see [20, 21] for a more complete discussion). For a given energy, we can record (over 
many different initial conditions) the coordinate and conjugate momentum of one particle 
every time the other particle has positive momentum and crosses a particular point. The 
resulting contours in phase space, collectively known as a Poincare section, display the 
transition from quasi-integrable to chaotic behavior in our system. 

For sufficiently low energies, the Poincare sections are given by two fixed points sur- 
rounded by a set of concentric contours. The fixed points correspond to motion in one of 
the two normal modes. It is useful to define the winding number w, which is the ratio of 
the number periods one particle completes for every full period completed by the other. At 
the linearized level away from the fixed point, the winding number w = u>i/u>2- If w is not a 
rational number, the trajectory will never quite return to its original position and thus fills 
one of the concentric contours. As the energy is increased, the winding number is detuned 
and eventually may even become rational. Hence, parts of the phase space acquire new fixed 
points with their own concentric contours. These correspond to nonlinear resonances. The 
last tori to break are those with the 'most' irrational w (the golden mean (V5- l)/2 is the 
'most' irrational number, as defined by the speed of convergence of its continued fraction 
expansion). The breaking of the original two islands into an increasing number is qualita- 
tively similar to the case of a double pendulum. Eventually, there is essentially no visible 
structure left in the Poincare section and we are in a regime of global chaos. We give an 



% = M^d^Hcd , 




1 





(6.3) 
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Figure 6.3: Poincare sections of collinear setup with K31 = 10, K32 = —10, K21 = —10 
and 6 3 = -1, 2 = 1 and energies E = {0.10,0.20,0.23,0.26,0.27,0.30}. Note that the 
k's form a closed loop. The horizontal axis represents the position of particle 3 while the 
vertical axis represents its conjugate momentum. Any given plot is produced by varying the 
initial positions and momenta of the two probes subject to a fixed total energy. The pair 
(x3(t),p3(t)) is plotted every time the resulting trajectory of particle 2 crosses some fiducial 
point {x2(t) = x c ) with positive momentum ip^if) > 0), i.e. roughly every time particle 2 
completes a full cycle as it oscillates back and forth. In the quasi-integrable regime, different 
initial conditions correspond to different contours. The first Poincare section shows a quasi- 
integrable behavior with two fixed points corresponding to the two low energy normal modes. 
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example of this in figure 6.3. Similar transitions to chaos are found for examples where the 
k's form closed and non-closed loops. In several examples where the k's form a closed loop 
and obey the triangle inequality, and the #'s have the same sign, the formation of islands 
around fixed points representing nonlinear resonances seems to be far less manifest in the 
Poincare sections. In other words, global chaos seems to set in much more quickly. We hope 
to study these issues systematically in the future. 

7 Trapping 

In this section, we envision a trapping problem. The setup consists of a localized bound state 
and another particle, which we take to be a probe, beginning inside the molecule. The probe 
begins its life at a random position well within the molecule. We explore the dynamical 
evolution of the probe as we vary the initial energy. 

7.1 Setup and Energetics 

Our setup will consist of a probe with charge 7 P = (1, 0) in the presence of a bulk molecule 
comprised of a number N c of fixed electric centers of charge 7 C = (0, k). The positions of the 
electric centers will be obtained by drawing random points from a ball of radius R mo i using 
the algorithm in [41, 42]. The classical probe Hamiltonian in this background is given by: 




Since all the background centers have the same charge, the all equal. Also, to 

ensure that trapping occurs we require that k and 9 P have opposite signs. The zero energy 
configurations are given by setting the second term in H pro b e to zero. As usual, there is a 
classical moduli space M. due to the fact that we have three probe coordinate degrees of 
freedom and we are solving only one equation. We could search for non-zero static minima 
of H pro be, but a simple computation of the gradient of the potential shows that there are 
only zero energy minima. The minimal energy required for the probe to reach infinity is 
Emin = 0p/2m p . 

For probe energies E in > 8 p /2m p the probe can easily escape the molecule. For energies 
in the range E in < 8p/2m p , we observe trapping. Our goal is to begin quantifying the amount 
of classical trapping. We do this by studying the fractional volume fv(E in ,t) covered by 
the probe as a function of initial energy E in and total trajectory time t. We estimate fy 
by studying how many centers the probe trajectory approaches to within one-half of the 
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average inter-particle distance r i p . ~ R mo i/N c ■ We again stress that we keep the molecule 
and initial positions of our probe fixed through all the trials, only varying the initial velocity 
of the probe. We take m p — 1, 9 P — —10, R mo i = 20, N c = 100 and k = 1 for the presented 
data. 

7.2 A trap 

At low energies, we witness characteristic trapping: figure 7.1 shows one such example 
where a probe is confined to less than 20% of the molecule, exploring the same part of the 
molecule over and over again. We remind the reader that it is possible to be trapped in one 
region indefinitely and this behavior should not (necessarily) be looked at as a failure of not 
integrating for a long enough period of time. Indeed, as is seen in the Euler-Jacobi flower of 
figure 5.1, probes can remain in one part of a molecule for arbitrarily long periods of time. As 
we increase the energy, we see a transition that opens up more of phase space to the probe. 
Little pockets in the potential landscape form through which the probe particle can escape 
and begin exploring other regions of the molecule. Often this happens by sudden jumps, as 
illustrated in the middle row of figure 7.1. We have tracked the energies of the probe and the 
numerics are stable. The jump is not due to an erroneous kick in the integrator but rather 
appears to be due to small pockets through which the probe can escape given enough time. 
Finally, at high energies (e.g. around one half of the escape energy), the probe uniformly 
explores the entire molecule, as illustrated in the last row of figure 7.1. 

It is also interesting to study the Schrodinger equation for the probe in this background 
to see if the wavefunction exhibits trapping via some avatar of Anderson localization. 4 

7.3 Topology of the potential landscape 

To illustrate the potential landscape and gain some intuition for the motion of trajectories, 
we set up a co-planar molecule. Again, all centers in the molecule have equal charge and their 
positions are chosen by uniformly selecting N c points on a disk of size R m oi [43]. A probe 
in this background will not remain in the plane due to the magnetic fields which will push 
it out. However, for a molecule where every center attracts the probe, at low energies the 
deviations from the plane are small relative to the size of the molecule, which can be made 
arbitrarily large. Thus, plotting equipotential contours over this two-dimensional molecule 
gives an accurate picture of the potential landscape which can be used to understand the 
trajectories. See figure 7.2 for such a comparison. For a fixed molecule size, as the magnitude 
of n is increased relative to the magnitude of 9, the topology of the equipotentials changes 

4 Thanks to Douglas Stanford for discussions on the quantum dynamics. 
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5000 10000 15 000 1 




Figure 7.1: The first row represents a low energy probe, which remains stuck in a subset of the 
phase space for seemingly arbitrarily long times. The second row represents an intermediate 
energy probe which illustrates the non-uniform escapes that occur from the low energy 
trapping behavior. We see that for a while it remains trapped in some subset of phase space, 
after which it escapes and gets stuck in some other subset of phase space. The final row 
represents a high energy probe which uniformly explores the molecule. The associated plots 
represent the percentage of the molecule explored as a function of the integration time, up to 
15000 time steps in increments of 1500. The tapering off of the high energy probe is simply 
due to saturating the entire molecule. Below these points the increase is very uniform. The 
initial energy increases from 50% of the escape energy in the first row to 60% of the escape 
energy in the third row. These percentages, however, are very dependent on the parameters 
(e.g. k, 9, etc.) in the problem. 
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Figure 7.2: These contour plots show equipotential surfaces in the plane of a 2D molecule 
consisting of one hundred centers. From left to right, we have chosen 

and in all cases 9 = — 10. We observe that as the magnitude of k increases, the minima (blue 
region), which initially lied near each center, are collectively expelled, forming an overall 
minimum that surrounds the molecule as a whole. For k — 1, the trajectory remains close to 
the plane of the molecule and has been superimposed on the left contour plot (transparent 
white line). The axes label the x and y positions of the probe particle. 

by expelling the low energy part of the landscape to the outside of the molecule, as can be 
seen in figure 7.2. This mimics the change in topology of the moduli space in going from 
Region V to Region IV in figures 4.1 and 4.2. Topology changes in the moduli space of the 
three center setup also occurs as k is increased while keeping 9 fixed as well as the distance 
between the background centers. 

8 Holography of Chaotic Trajectories? 

We end our journey by discussing how the picture we are developing may fit into the broader 
context of holography. The usual interpretation of a large black hole in an asymptotically 
anti-de Sitter space is that we have prepared the dual CFT in some finite temperature 
state. On the other hand, the presence of a vast number of distinct entropically relevant 
multicentered black hole configurations inside an anti-de Sitter universe [45] implies a vast 
number of minima in the free energy of the dual CFT (as a function of configuration space). 
In particular, the usual assumptions of the no-hair theorem fail since a set of macroscopic 
charges does not uniquely fix the bulk solutions. In fact, the solutions are characterized 
by a large collection of multipole moments. Furthermore, though most of them do not 
constitute true ground states, they can be very long lived, decaying mostly through thermal 
and quantum tunneling. 

We can access information about the relaxation and response of the CFT by computing 
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boundary-to-boundary correlators in the bulk. In the large frequency (or large mass) limit 
two-point functions have been associated to bulk geodesies which begin and end their life 
near the boundary of AdS [46, 47, 48, 49]. Such geodesies will become highly complex and 
chaotic in the bulk due to the presence of the non-trivial black molecule, as evidenced by 
our simpler setup. In fact, a geodesic may become trapped in some very long lived unstable 
orbit before escaping back to the boundary. Thus, the two-point function expressed as a path 
integral over bulk trajectories and the applicability of the saddle point approximation may 
be a somewhat involved issue. This picture suggests that the linear response properties of 
the dual CFT, to the extent that they are captured by the two-point function in the geodesic 
approximation, in the multicentered/glassy phase are rather different from those in a usual 
thermal state, where for example the motion of geodesies is integrable. The motion of a very 
massive probe or high energy graviton falling into the bulk corresponds to a point-like source 
cascading to lower energies (and covering larger size) in the CFT. Eventually the excitation 
returns back to a point-like source at some other point on the two-sphere where the CFT 
resides. From the bulk point of view this is when the particle dropped into the molecule 
returns back out. The chaotic nature of the bulk physics suggests chaotic behavior of the 
boundary theory itself. 

One may also consider the dynamical features at zero temperature for which asymptot- 
ically AdS3 multicentered configurations are known [50]. The possible presence of classical 
chaotic behavior of the bulk AdS3 should correspond to quantum dynamics in the dual CFT2. 
One effect of particular interest in chaotic systems is known as quantum scarring, where it has 
been observed that the wavefunction of a chaotic system peaks on closed classical trajectories 
[51]. We hope to explore these issues further in future work. 
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A Two-Body Problem 



In this appendix we discuss the motion of a single probe particle in the background of a 
fixed charge sitting at the origin. This was studied at length in [19]. The Hamiltonian of 
this system is given by 

H= L<»- A -r + L& +t Y • P = mx+A (A1) 

and is conserved. For simplicity we choose k > and allow 9 to be either positive or 
negative. Other than the Hamiltonian, this system admits two vector-valued conserved 
quantities known as the angular momentum L and the Runge-Lenz vector n. Explicitly 

L = x x (p - A) + — x , and n = f x + - L J x (p - A) . (A.2) 
2r \ 9 J 

Since this system is superintegrable, the probe particle's trajectories can be found alge- 
braically. 

First, notice that n-x = implying that n is perpendicular to the plane of motion of the 
probe particle. We use this fact to orient our axes such that n = |n|5. It is straightforward 
to show that 

which implies that |L| > k/2. With this choice of coordinates, the particle's trajectory is 
constrained to lie in a plane of constant z. The magnitude of z can be obtained by computing 
n ■ x = \n\z = — (h 2 — ^ /9, giving 



We have yet to choose an orientation for the x — y plane; we do so by aligning our coordinates 
such that L y = and L x points in the positive x direction. The components of the angular 
momentum are given by 



L x = x \\?- — (l*--), L z = J—(v--). (A.5) 
1 2mH V 4 J ' y 2mH \ 4 J K J 

We can determine the particle's trajectory explicitly by noticing that (A.2) implies that 
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Figure A.l: Scatering angle. 



L-x = Kr/2or 

(1 - e 2 ) (x - x ) 2 + 2ie[x - x ) + y 2 = t 2 , (A.6) 

which is the equation for a conic section in cartesian coordinates. The quantities e and I are 
the eccentricity and the semi-latus rectum of the conic section respectively and are given by 




2mH V 4 J ' nV8mH ' 



The quantity x = y L 2 — (L 2 — ^)/(|^| _ \/2mH) is the location of one of the foci of 
the conic section. 

The overall shape of a conic section is determined by its eccentricity, with elliptic orbits 
corresponding to e < 1, while parabolic and hyperbolic orbits correspond to e = 1 and e > 1 
respectively. Intuition predicates that bound orbits should only happen for 9 < 0, while 
6 > gives rise to parabolic or hyperbolic orbits. For positive 6 the Hamiltonian is bound 
such that 2mH > 6 2 which implies that e > 1, thus verifying our intuition. 

For e < 1 the length of the semi-major axis a of the elliptic orbit is given by 

2t nJH/2m ^ . . . 

V +0(1) (A.8) 



1 e -^escape H 



so as the energy approaches the escape energy H escape = 9 2 /2m (or as e approaches 1), the 
size of the bound orbit diverges. We end this appendix with a discussion of scattering for 
e > 1. Since the trajectory of the particle is given by a conic section in the x — y plane, the 
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scattering angle as defined in figure A.l is given by 



(p = 2 arccot (Ve 2 - 1 j . (A.9) 

For parabolic orbits (e = 1), (p — tt and we see that the particle completely back scatters. 
As we increase e the scattering angle decreases monotonically. 
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